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This talk brings attention to what is knowable from perturbative QCD theory on 
two-parton distribution functions in the light of CDF mcasuruments of the inclusive 
cross section for double parton scattering. 

The Collider Detector at Fermilab (CDF) Collaboration has measured a large number 
of double parton scattering [Tj providing new and complementary information on the 
structure of the proton and parton-parton correlations. The possibility of observing 
two separate hard collisions has been proposed since long 2 , and from that has also 
developed in a number of works [31 0] . A brief review of the current situation and some 
progress in the modeling account of correlated flavour, colour, longitudinal and transverse 
momentum distributions can be found in Ref . [5] . Multiple interactions require an ansatz 
for the structure of the incoming beams, i.e. correlations between the constituent partons. 
As a simple ansatz, usually, the two-parton distributions are supposed to be the product 
of two single-parton distributions times a momentum conserving phase space factor. In 
recent papers [B] it has been shown that this hypothesis is in some contradiction with 
the leading logarithm approximation of perturbative QCD (in the framework of which a 
parton model, as a matter of fact, was established in the quantum field theories [7]). 

In order to be clear and to introduce the denotations let us recall that, for instance, 
the differential cross section for the four-jet process (due to the simultaneous interaction 
of two parton pairs) is given by [5] 

da = y^J^^^^^—^ Dp{xi,X3) Dp{x2,X4), (1) 

, CTeff 
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where ddy stands for the two-jet cross section. The dimensional factor <7cS in the denom- 
inator represents the total inelastic cross section which is an estimate of the size of the 
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hadron, (Toff — 27rrp (the factor 2 is introduced due to the identity of the two parton pro- 
cesses). With the effective cross section measured by CDF, (o'off)cDF = (14.5 ± 1-7^2 3) 
mb [1], one can estimate the transverse size ~ 0.5 fm, which is too small in comparison 
with the proton radius Rp extracted from ep elastic scattering experiments. The rela- 
tively small value of ((Toff)cDF with respect to the naive expectation 27ri?p was, in fact, 
considered as evidence of nontrivial correlation effects in transverse space. But, apart 
from these correlations, the longitudinal momentum correlations can also exist and they 
were investigated in Ref. [5]. The factorization ansatz is just applied to the two-parton 
distributions incoming in Eq. (fT|): 



Dp{xi,Xj) = Dp{xi,Q'^) Dp{xj,Q'^) (1 



(2) 



where Dp{xi,Q^) are the single quark/gluon momentum distributions at the scale 
(determined by a hard process) . 

However many parton distribution functions satisfy the generalized Gribov-Lipatov- 
Altarelli-Parisi-Dokshitzcr (GLAPD) evolution equations derived for the first time in 
Refs [SI [9] as well as single parton distributions satisfy more known and cited GLAPD 
equations [7l [10] . Under certain initial conditions these generalized equations lead to so- 
lutions, which are identical with the jet calculus rules proposed originally for multiparton 
fragmentation functions by Konishi-Ukawa-Veneziano |llj and are in some contradiction 
with the factorization hypothesis ^ . Here one should note that at the parton level this 
is the strict assertion within the leading logarithm approximation. 

After introducing the natural dimensionless variable 
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where g{fJ-^) is the running coupling constant at the reference scale /i^, Uf is the number 
of active flavours, Aqcd is the dimensional QCD parameter, the GLAPD equations 
read [3 [10] 



dDl{x,t) 
dt 



(3) 



They describe the scaling violation of the parton distributions Dj{x,t) inside a dressed 
quark or gluon (i, j = q/g). 
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We will not write the kernels P explicitly and derive the generalized equations for 
two-parton distributions Dj^''^(xi,X2,t), representing the probability that in a dressed 
constituent i one finds two bare partons of types ji and j2 with the given longitudinal 
momentum fractions xi and X2 (referring to [H [71 [H [9l [TOj for details), we note only 
that their solutions can be represented as the convolution of single distributions [HI [S] ■ 
This convolution coincides with the jet calculus rules as mentioned above and is 
the generalization of the well-known Gribov-Lipatov relation installed for single func- 
tions [3 (the distribution of bare partons inside a dressed constituent is identical to the 
distribution of dressed constituents in the fragmentation of a bare parton in the leading 
logarithm approximation). The obtained solution shows also that the double distribution 
of partons is correlated in the leading logarithm approximation: 

Di^^-^{x^,X2,t) ^ Di^{xut)Df{x2,t). (4) 

Of course, it is interesting to find out the phcnomenological issue of this parton level 
consideration. This can be done within the well-known factorization of soft and hard 
stages (physics of short and long distances). As a result the equations ([3]) describe the 
evolution of parton distributions in a hadron with t (Q^), if one replaces the index i by 
index h only. However, the initial conditions for new equations a,t t — (Q^ = are 
unknown a priori and must be introduced phenomenologically or must be extracted from 
experiments or some models dealing with physics of long distances [at the parton level: 
Dl{x,t = 0) = SijS{x-l); D{'^^xi,X2,t = 0) = 0] . Nevertheless the solution of the 
generalized GLAPD evolution equations with the given initial condition may be written 
as before via the convolution of single distributions [H [9] . This result shows that if the 
two-parton distributions are factorized at some scale then the evolution violates this 
factorization inevitably at any different scale {Q^ ^ /i^), apart from the violation due to 
the kinematic correlations induced by the momentum conservation. 

For a practical employment it is interesting to know the degree of this violation. 
Partly this problem was investigated theoretically in Refs. [SI [12] and for the two-particle 
correlations of fragmentation functions in Ref. [13| . That technique is based on the Mellin 
transformation of distribution functions and the asymptotic behaviour can be estimated. 
Namely, with the growth of t (Q^) the correlation term becomes dominant for finite xi 
and X2 [E] and thus the two-parton distribution functions "forget" the initial conditions 
unknown a priori and the correlations perturbatively calculated appear. 

The asymptotic prediction "teaches" us a tendency only and tells nothing about the 
values of xi, X2,t{Q'^) beginning from which the correlations are significant. Naturally 



A.M. Snigirev 



HSQCD 2008 



numerical estimations can give an answer to this specific question. We do it using the 
CTEQ fit [14] for single distributions as an input. The nonperturbative initial conditions 
Dj^{x,0) are specified in a parametrized form at a fixed low-energy scale Qq = ii = 
1.3 GeV. The particular function forms and the value of Qq are not crucial for the 
CTEQ global analysis at the flexible enough parametrization. The results of numerical 
calculations are presented in Fig. 1 for the ratio: 



Rix,t) = (^DllQ^^^{xuX2,t)/DS{xi,t)Da{x2,t){l-x,^X2)^) 



X1—X2—X 



(5) 
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Figure 1: The ratio of perturbative QCD correlations to the factorized component for 
the double gluon-gluon distribution in the proton as a function of x = xi — X2 for three 
values of Q = 5 (solid), 100 (dashed), 250 (dash-dotted) GeV. 



Figure 1 shows that at the scale of CDF hard process (^ 5 GeV) the ratio ([5]) 
is nearly 10% and increases right up to 30% at the LHC scale 100 GeV) for the 
longitudinal momentum fractions a; < 0.1 accessible to these measurements. For the 
finite longitudinal momentum fractions x ~ 0.2 0.4 the correlations are large right up 
to 90% . They become important in more and more x region with the growth of t in 
accordance with the predicted QCD asymptotic behaviour [T^. 
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The correlation effect is strengtliened insignificantly (up to 2%) for the longitudinal 
momentum fractions x < 0.1 when starting from the slightly lower value Qq ~ 1 GeV 
(early used by CTEQ Collaboration). We conclude also that R{x,t) const at a; ^ 
most likely, calculating this ratio (~ 0.1) at a^min — 10^^. 

Seemingly the correction to the double gluon-gluon distributions at the CDF scale 
can be smoothly absorbed by uncertainties in the (leff increasing the transverse effective 
size rp by a such way. But this augmentation is still not enough to solve a problem 
of the relatively small value of rp with respect to the proton radius without nontrivial 
correlation effects in transverse space W . 

Recently a nonminor role of the QCD evolution of multiparton distribution functions 
has been also demonstrated [Hj. In the case of multiple production of W bosons with 
equal sign, the terms with correlations may represent a correction of the order of 40% of 
the cross sections, for pp collisions at 1 TeV cm. energy, and a correction of the order 
of 20% at 14 TeV. In the case of bb pairs the correction terms are of the order of 10-15% 
at 1 TeV and of the order of 5% at 14 TeV. 

In summary, the numerical estimations show that the leading logarithm perturbative 
QCD correlations are quite comparable with the factorized distributions. With increasing 
a number of observable multiple collisions (statistic) the more precise calculations of their 
cross section (beyond the factorization hypothesis) will be needed also. In order to obtain 
the more delicate their characteristics (distributions over various kinematic variables) it 
is desirable to implement the QCD evolution of two-parton distribution functions in some 
Monte Carlo event generator as this was done for single distributions. 
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